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Resumen 


Balocchi, F., Pizarro, R., Morales, C., € Olivares, C. 
(septiembre-octubre, 2014). Modelamiento matemático de 
caudales recesivos en la región mediterránea andina del 
Maule; el caso del estero Upeo, Chile. Tecnología y Ciencias del 
Agua, 5(5), 179-188. 


En el presente estudio se analizaron cuatro modelos 
matemáticos de estimación de caudales recesivos, con 
el fin de determinar cuál de ellos obtiene los mejores 
resultados en una cuenca andina del Maule, zona 
mediterránea de Chile central. Esto se realizó a través 
del análisis de 25 crecidas en la época estival 1971-2003, 
en donde se consideró el aporte exclusivo de las aguas 
subterráneas a partir del tercer punto de quiebre de 
la curva de bajada del hidrograma. Los resultados 
obtenidos demuestran que el modelo exponencial 3 es 
el que presenta las mejores estimaciones. En el caso del 
mejor lapso para el ajuste de los modelos, el tiempo de 
240 horas es el que obtiene los mejores resultados, según 
el test de Bland y Altman. 


Palabras clave: hidrograma, caudales recesivos, curva de 
bajada, modelos matemáticos. 


Introducción 


La predicción de los caudales recesivos es 
una metodología ampliamente usada en la 
planificación y en la investigación hidrológica 
(Smakhtin, 2001), por lo que se transforma 
en una importante herramienta hidrológica, 
especialmente en periodos estivales, ya que 
en esta época la población es alimentada por 
las reservas de la cuenca (Wittenberg, 1999). 
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This study analyzed four mathematical models for recessive flow 
estimation in order to determine which one would provide the best 
results for an Andean basin in Maule, central Chile. This was 
accomplished through the analysis of 25 summer floods occurring 
between 1971 and 2003, which identified groundwater as the 
exclusive supply as of the third inflection point of the falling limb 
of the storm's hydrograph. The results indicate that the exponential 
model 3 provides the best estimates. According to the Bland and 
Altman agreement test, the best model adjustment period was 240 
hours. 
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Esta situación deriva en que el conocimiento 
adquirido sobre las reservas de agua de una 
cuenca es de vital importancia para el desarrollo 
humano y productivo de un país. 

Es por esto que surge la necesidad de 
estimar la disponibilidad de este recurso para 
su utilización en los tiempos estivales. Así, 
la modelación matemática constituye una 
importante herramienta para el análisis y una 
mejor comprensión del comportamiento del 
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agua en estos periodos, considerando que tales 
modelos pueden estudiar el comportamiento 
de sistemas complejos ante situaciones difíciles 
de observar en la realidad (Sujono, Shikasho, éz 
Hiramatsu, 2004). Estrela (1992) señala que los 
modelos matemáticos para estimar caudales 
constituyen una forma válida para resolver los 
problemas hídrico-teóricos y, además, son los 
que tienen un mayor sentido físico-hidrológico. 

Por otra parte, Linsley, Kohler y Paulhus 
(1949), y Pizarro (1993) señalan que si se conoce 
la forma de la curva de recesión del hidrograma 
es factible calcular el flujo de corrientes durante 
lapsos de sequía o estiaje. Así, el hidrograma 
muestra el comportamiento de los caudales 
con respecto al tiempo, el cual está compuesto 
por una curva de concentración, una curva 
de descenso y una curva recesiva (Maidment, 
1992; Pizarro, 1993; Aparicio, 2003; Bedient 
Huber, 2002; Brodie, Hostetler, €z Slatter, 2007). 
La curva recesiva es aquella donde el flujo 
proviene exclusivamente de aguas subterráneas 
y permite calcular la disponibilidad de agua 
en los periodos entre crecidas, lo cual es de 
particular importancia en los lapsos estivales. 
En la figura 1 se puede observar que la curva 
recesiva O de agotamiento comprende desde el 
punto C en adelante (Linsley et al., 1949); para 
el caso del segundo punto de quiebre y el punto 
C', indica el inicio de los caudales recesivos, 
según la metodología del tercer punto de 
quiebre (Pizarro et al., 2013), esta última 
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Figura 1. Hidrograma de crecida. 


utilizada en el presente estudio. Lo anterior se 
debe a que al asumir el tercer punto de quiebre 
como inicio de la curva recesiva se adquiere 
una mejor estabilidad de los datos de caudal, 
reduciendo la variabilidad de los mismos, lo 
cual permite mejores niveles de ajuste. 

La separación del hidrograma propuesta 
por Linsley et al. (1949) y corroborada por 
Pizarro (1993) es la más utilizada en Chile. 
Dicha separación consiste en dibujar la curva 
de bajada del hidrograma en un gráfico 
semilogarítmico, en donde el segundo punto 
de quiebre corresponde al punto de inicio del 
hidrograma, donde sólo existen aportes de 
los caudales recesivos (Vilaró, 1976; Bedient éz 
Huber, 2002). Existe otro tipo de metodologías 
para la separación del hidrograma, como la de 
Bedient y Huber (2002), y la de Viessman y Lewis 
(2003), ofreciendo buenas aproximaciones, 
que en general no presentan variaciones de 
relevancia al momento de determinar el punto 
de inicio de los caudales recesivos (Ponce, 
1989). En el ámbito académico, es la teoría 
de Hewlet y Hibbert (1967) la más aceptada 
respecto a la separación de hidrogramas. Esta 
metodología señala que el agua nueva de lluvia 
empuja el agua antigua del suelo y esta última 
sería la que se une al caudal superficial y/o 
subterráneo. Sin embargo, Brooks, Barnard, 
Coulombe y Mcdonnell (2009) demostraron 
que no necesariamente se cumple dicho 
comportamiento. Por último y para determinar 
el comportamiento de los diferentes modelos 
matemáticos recesivos, se optó por mantener 
de manera constante la metodología de 
separación de Linsley. 

Considerando lo anterior y con base en 
el conocimiento adquirido de las curvas 
recesivas, se modeló el comportamiento en 
periodos estivales de los caudales subterráneos 
que afloran en el estero Upeo, río Lontué, Chile, 
con el fin de encontrar el modelo que mejor se 
ajuste a la zona del estudio. Esto, debido a la 
disminución de la disponibilidad hídrica en 
la zona mediterránea de Chile, en especial en 
época estival, lo que sumado a la determinación 
del modelo correcto en la cuenca en estudio, 
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garantizaría una adecuada estimación de los 
caudales recesivos y el posterior cálculo de las 
reservas de la cuenca. Así, la estación utilizada 
en el presente estudio, estero Upeo en Upeo, 
está ubicada en la coordenada UTM 309 573 E— 
6 105 697 N (WGS84/h19). 

La zona de estudio (figura 2) posee una 
precipitación de 1 800 mm/año y corresponde 
a la subcuenca del río Lontué, ubicada entre 
309 340 E — 6 095 942 N y 342 916 E- 6 114 754 
N (WGS84/h19), abarcando una superficie de 
2 784 km? que, en conjunto con el río Teno, 
ambos de régimen nivo pluvial y originados en 
la cordillera de los Andes, dan como origen al 
río Mataquito. El río Lontué recibe por su parte 
norte el aporte de dos esteros tributarios: Upeo 
y Chenquelmo. 

Las temperaturas registradas en dicha 
estación fluctúan entre los 20 y 4 “C. En cuanto 
a sus caudales, el máximo registrado en la 
estación corresponde a 490 m*/s en el año 2000. 

El rasgo geomorfológico principal de 
la zona corresponde a la cordillera de los 
Andes. En cuanto a la geología, el tipo de roca 


predominante es el paleozoico mesozoico y en 
algunos sectores el tipo cenozoico. Además, 
el sector se caracteriza por la presencia de 
sedimentos glaciares del tipo morrena (CNR, 
1978). 


Metodología 


Para la selección de la estación, se tomó como 
base que ésta presentara un régimen natural 
de caudal, monitoreado por una estación 
fluviométrica que presentara un limnígrafo. La 
estación seleccionada posee datos desde 1964 
hasta la fecha. Por tanto, se escogió la estación 
Estero Upeo, en Upeo, cuenca del Mataquito, 
región del Maule (figura 2). 

Así, la información requerida consideró 
las curvas de descarga y los limnigramas 
desde 1964 hasta 2003, correspondientes a 
la Estación Estero Upeo, en Upeo, la cual fue 
proporcionada por la Dirección General de 
Aguas (DGA), institución oficial encargada 
de la administración y medición del recurso 
hídrico en Chile. 
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Figura 2. Ubicación de la zona de estudio. 
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Con la información recolectada se selec- 
cionó el periodo que poseía crecidas que 
fuesen claramente identificables en época 
estival (octubre-abril). De acuerdo con esto, 
se seleccionaron 25 crecidas entre los años 
hidrológicos 1971-2003, con las que se procedió 
a construir el hidrograma de crecida de cada 
tormenta. Cabe recordar que se seleccionaron 
las últimas crecidas de cada año, casi al final 
del año hidrológico, con el fin de asegurar 
la presencia de caudales recesivos en lapsos 
importantes. 

Los lapsos mencionados corresponden 
al tiempo que transcurre desde el inicio de 
un caudal recesivo (Q,) y hasta el inicio de 
una nueva crecida. Además, estos lapsos son 
esenciales para el cálculo de los parámetros 
de los modelos utilizados. En este contexto, el 
menor periodo de crecida encontrado fue de 11 
días y el máximo de 239 días. 

En este contexto, el presente estudio 
consideró cinco lapsos correspondientes a 48 
horas, 72 horas, 168 horas, 240 horas y el total 
de horas de cada crecida seleccionada, debido 
a que existe la hipótesis de que a mayor tiempo 
transcurrido desde la coordenada inicial (Q,, t.), 
el ajuste de cada modelo debería ser de mejor 
calidad. Así, las coordenadas para el lapso t = 
48 horas sería (Q,, t.), (Q yyy t,5); para t= 72 horas 
sería (Qy t), (Q,, t,,), y así sucesivamente, 
para cada lapso escogido. No obstante, y por 
su representatividad en el estudio, sólo se 
consideraron los lapsos de 240 y el total de 
horas de cada crecida. 

Por otra parte, los modelos matemáticos 
empleados correspondieron a cuatro tipos: uno 
potencial y tres exponenciales, en los cuales los 
parámetros utilizados son el caudal inicial (Q ,), 
el coeficiente de agotamiento (a), un parámetro 
particular (n) del modelo exponencial 2, la 
constante de Neper (e) y los tiempos t y t, 
donde t, corresponde al tiempo inicial de los 
modelos. 
modelos se basa 
considerando el valor de caudal inicial de la 


El ajuste de estos 


recesión (Q,), que se produce al tiempo t,= 0; 
ello involucra considerar uno o más puntos en 


el instante t, + dt, que define otra coordenada Q, 
con lo cual es posible obtener los parámetros de 
la ecuación. Por ello es importante estudiar si 
el diferencial de tiempo incide en la calidad de 
predicción de los caudales recesivos. Con esto, 
los modelos empleados fueron: 


Modelo potencial: Q(t) = Q,( + a : t)? 
(Cirugeda, 1985) (1) 


Modelo exponencial 1: Q(t) = Q, : exo 
(Remenieras, 1971) (2) 


Modelo exponencial 2: Q(t) = Q, : ect” 
(Cirugeda,1985) (3) 


Modelo exponencial 3: Q(t) = Q, : e 
(elaboración propia) (4) 


Para el caso del modelo exponencial 3, 
el cual se propone como un nuevo modelo 
para la zona en estudio, se generó de una 
serie de modelos construidos y simulados. Su 
determinación matemática fue la siguiente: 


Se sabe que de < 0 donde Q = caudal y t = tiempo 


Si se asume que la variación de caudal en 
el tiempo es proporcional al mismo caudal e 
inversamente proporcional a la raíz del tiempo, 
se tiene: 


a = -a( E operando S e -a( 5, 


Se sabe además que cuando Q = Q,, se define 
que t=0 = £,. Integrando entre los límites t y £, 
se tiene: 

In QS" =aWt $, lo cual queda: 

InQ-InQ, = -2a( ye - Vio). Por ende: 
2h 

Q(+) =0Q, sl 1) 


Para determinar los parámetros de los 
modelos y así poder ajustar éstos a la curva 
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recesiva, el procedimiento implicó considerar 
una primera coordenada, que es la del inicio de 
los caudales recesivos (t, Q,). Posteriormente 
se consideró una segunda coordenada, la cual 
permite calcular el coeficiente de agotamiento 
de los modelos. Sin embargo, el modelo 
exponencial 2 requiere considerar una tercera 
coordenada para satisfacer el sistema de 
ecuaciones. 

La validación de los modelos se realizó 
a través de distintas pruebas e indicadores 
estadísticos no paramétricos. Esta etapa de 
validación estadística consideró los resultados 
de las 25 crecidas seleccionadas, para los cuatro 
modelos propuestos. En el marco descrito se 
utilizaron los siguientes indicadores y pruebas 
estadísticas: 


a) Coeficiente de determinación (R?) 


Este coeficiente, también denominado Nash- 
Sutcliffe, determina qué proporción de la 
variación total de los caudales observados 
(variable independiente), es explicada por los 
caudales modelados (variable dependiente). La 
forma de determinarla es la siguiente: 


yi: caudales observados 


AZ 
R?=1- Ali) yi: promedio de los 
Ss caudales observados 


Yi: caudales modelados (5) 


b) Error estándar de estimación (EEE) 


Determina la disparidad promedio entre los 
caudales observados y los estimados, donde 
los valores cercanos a O indican que el modelo 
hace una buena descripción de los caudales de 
acuerdo con la fórmula del indicador: 


y: caudales observados 
y: caudales modelados 


n: número de datos de 
caudal (6) 


c) Prueba U de Mann-Whitney 


Esta prueba determina si dos muestras 


independientes provienen de la misma 
población (Montgomery $ Runger, 1996; 
Mendenhall € Sincich, 1997). Asimismo, al 
ser una prueba no paramétrica, no requiere 
de los supuestos de normalidad e igualdad de 
varianzas (Mason é Lind, 1995). 

Se establecen la hipótesis nula y la hipótesis 


alternativa: 


H,: las distribuciones de las variables R, y R, 
son idénticas. 

H,; las distribuciones de las variables R, y R, 
no son idénticas. 


Donde R, corresponde a los datos reales y R, 
a los datos modelados: 


n,-n,+1 ends 
Ua=(n,-1,)+ (Pza - SR, estadístico 
para muestras pequeñas (1, y n, < 25) (7) 


DE SR, 8 -|(n, e) (Bat) 


fuen) 


dístico para muestras grandes (n, y n,>25) (8) 


esta- 


Donde: 

n,: tamaño de la muestra RR; n,: tamaño de 
la muestra R,; Y R,: suma de los rangos para la 
muestra R,; 2 R,: suma de los rangos para la 
muestra R,. 

Regla de decisión: se rechaza H,, a un nivel 
de significancia al, si: 


Z> Za] (9) 


Donde Za: valor crítico obtenido de la tabla 
normal estándar, a un nivel de significancia de 
a del orden de 0.05 y 0.01. 
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d) Test de concordancia de Bland y Altman 
(ACBA) 


Este test determina si los modelos de estimación 
concuerdan lo suficiente como para que puedan 
declararse intercambiables (Carrasco ézx Jover, 
2004; Cox, 2006). En términos estadísticos, se 
procede a calcular las diferencias promedio (dp) 
entre el valor real y el estimado, cuya diferencia 
representa el error sistemático, mientras que 
la varianza de las diferencias (DS) mide la 
imprecisión. En este marco se establecen los 
límites de concordancia del 95% (Carrasco 
éx Jover, 2004; Paradis, Lefebvre, Morin, « 
Gloaguen, 2010). Este límite está definido por 
la expresión: 


LC =dp=2-DS (10) 


Así, la relación valor real versus valor 
estimado con un dp en valor absoluto lo más 


cercano a cero será la de mejor ajuste. Si llegase 
a existir una igualdad o una mínima diferencia 
del valor dp, el mejor ajuste corresponderá al 
que presente un menor valor DS y sus límites 
de concordancia serán más estrechos (Bland é- 
Altman, 1999). 


Resultados y discusión 


A modo de ejemplo visual, en la figura 3 se 
presentan dos crecidas diferentes con sus 
respectivos modelos (ajustado para el total de 
datos) y caudales reales, con el fin de mostrar 
la alta variabilidad de los modelos. 

En el cuadro 1 se muestran los coeficientes 
de agotamiento (a) promedio para los cuatro 
modelos propuestos para los cinco lapsos 
establecidos y para las 25 crecidas. 

Del cuadro 1 se desprende que el valor del 
coeficiente de agotamiento va disminuyendo a 
medida que aumenta el tiempo de estimación 
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Crecida 21/12/1971 
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Figura 3. Ejemplo gráfico de dos crecidas con los cuatro modelos utilizados. 
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Cuadro 1. Coeficiente de agotamiento (a). 


a (240 h) a total (*) 
Promedio 0.0021 0.0012 
Modelo potencial Coef. variación 0.7437 0.8070 
Promedio 0.0031 0.0015 
Modelo exponencial 1 Coef. variación 0.6281 0.9102 
Promedio 0.0244 0.0224 
Modelo exponencial 2 Coef. variación 0.6281 0.4604 
Promedio 0.0244 0.0224 
Modelo exponencial 3 Coef. variación 0.6269 0.4601 


(*) Considera la totalidad de datos con que se cuenta; a, coeficiente de agotamiento. 


Cuadro 2. Resumen resultados test estadísticos R?, EEE y Mann-Whitney, aplicados a cada modelo. 


Test R? EEE U Mann-Whitney 
Modelo a (240 h) a total (*) a (240 h) a total (*) a (240 h) a total (*) 
Potencial 0:22 0.49 1.27 1.45 4% 20% 
Exponencial 1 0.08 0.3 1.47 1.83 4% 20% 
Exponencial 2 0.16 0.53 1.51 1.24 0% 32% 
Exponencial 3 0.28 0.59 1.2 11 0% 44% 


(*) Considera la totalidad de datos con que se cuenta; a,, coeficiente de agotamiento. 


debido a que al aumentar el tiempo de ajuste 
(diferencial) disminuye la pendiente entre el 
primer y el último punto. 

Para el modelo exponencial 2, el coeficiente 
de agotamiento fue en aumento a medida que 
se incrementaron los tiempos de estimación. 
Esto se explica, ya que en su fórmula de cálculo, 


Mon” 


al estar elevado a “n” el tiempo t, amplifica el 
valor del mismo, siendo este modelo el que 
registró los valores más altos en sus coeficientes, 
en comparación con los otros tres modelos. 

En cuanto al modelo exponencial 3, 
demostró una estabilidad en los valores del 
coeficiente, además de que es el modelo con 
el menor coeficiente de variación, situación 
que podría demostrar ser un buen modelo de 
ajuste. 

Por otro lado y en relación con el parámetro 


Mon 


n” del modelo exponencial 2, los datos 
denotan un aumento a medida que se van 
incrementando los tiempos de estimación: 
0.4834 para 240 horas y 0.4940 para el total de 
horas de la crecida (valores promedio). En este 


sentido, los resultados obtenidos concuerdan 


con el comportamiento del parámetro «a. Esto 
puede deberse a que al aumentar el parámetro 


Mon 


n”, el coeficiente de agotamiento (a) también 
lo hace (por fórmula de cálculo). 

En el cuadro 2 se presenta el resumen de 
los resultados del coeficiente de determinación 
(R2) promedio, del error estándar de estimación 
promedio (EEE), del test U Mann-Whitney y del 
test de Bland-Altman para los cuatro modelos 
propuestos, para los dos lapsos considerados y 
para las 25 crecidas. 

El cuadro 2 muestra que los mayores valores 
promedio de R? se encuentran concentrados en 
el total de datos. Esta situación podría deberse a 
que la fórmula de cálculo del R?es muy sensible 
a los cambios en el denominador. No obstante, 
los valores del coeficiente de determinación 
fueron muy bajos y no parecen representar de 
manera adecuada el comportamiento de los 
modelos matemáticos utilizados. Asimismo, 
y para el total de datos con que se cuenta, el 
modelo exponencial 3 obtiene el coeficiente más 
alto, aunque los valores son bajos en general. 


Por otro lado, el modelo potencial obtuvo los 


Tecnología 
E 


Ciencias del Agua, vol. V, núm. 5, septiembre-octubre de 2014 


y 


gua, vol. V, núm. 5, septiembre-octubre de 2014 


Ciencias del A 


d vy80]0u92] 


Balocchi et al., Modelamiento matemático de caudales recesivos en la región mediterránea andina del Maule; el caso del estero Upeo, Chile 


mejores resultados de R”, seguido del modelo 
exponencial 3. 

Para los modelos exponencial 1 y 2 se 
visualizan los coeficientes promedios más 
bajos y se observó una irregularidad en sus 
valores, no encontrándose una tendencia 
clara. Así, el modelo exponencial 1 fue el que 
obtuvo los valores más bajos del coeficiente 
de determinación, lo que señalaría que este 
modelo no parece ser adecuado para su uso. 

Para el caso del EEE (cuadro 2), se 
puede observar que los menores valores 
promedio (el 80% de ellos) se concentran en 
el modelo exponencial 3. En cuanto al modelo 
exponencial 1, fue el que obtuvo los valores 
más altos del EEE, lo cual estaría demostrando 
que este modelo no representaría de forma 
fehaciente el comportamiento de los caudales 
recesivos de la cuenca en estudio. 

También se puede comentar que los valores 
del EEE disminuyen a medida que aumentan 
los tiempos de estimación, con lo cual podría 
confirmarse la hipótesis de que a mayor tiempo 
transcurrido, mejor es el ajuste del modelo. 

Los resultados obtenidos en el test U de 
Mann-Whitney, calculados con un nivel de 
significancia del 5% (cuadro 2), muestran que 
en el total de datos se concentra el mayor 
porcentaje de pruebas aprobadas, observándose 
que el modelo exponencial 3 es el que obtiene 
los mejores resultados, seguido muy de cerca 
por el modelo potencial. En cuanto al modelo 


exponencial 1, de nuevo se descartó como 
predictor de caudales recesivos, al obtener los 
resultados menos satisfactorios en este test. 

Por último, los bajos resultados obtenidos 
en la prueba U de Mann-Whitney pueden 
explicarse debido a que el nivel de significancia 
utilizado posee un área de aceptación pequeña, 
por lo que se convierte en una prueba muy 
estricta. 

En el caso del test de Bland y Altman 
(ACBA), cabe señalar que se escogió teniendo 
en consideración el hecho de que tanto sus 
creadores como otros investigadores discre- 
pan de la confiabilidad del coeficiente de 
determinación R? ya que si dos instrumentos 
miden sistemáticamente cantidades diferentes 
uno del otro, la correlación puede ser perfecta, 
a pesar de que la concordancia sea nula (Pita éz 
Pértegas, 1998). 

En el cuadro 3 se encuentran los resultados 
obtenidos para las diferencias promedio (dp) y 
para la varianza de las diferencias (DS) en los 
cuatro modelos utilizados. 

Según el cuadro 3 y observando las 
diferencias promedio, el mejor modelo es el 
exponencial 2, seguido del modelo exponencial 
3. Sin embargo, al observar la desviación 
estándar de las diferencias promedio, es posible 
detectar una clara supremacía del modelo 
exponencial 3. 

De nuevo se descarta como modelo 
predictor al modelo exponencial 1 debido a 


Cuadro 3. Diferencias promedio (dp) y varianza de las diferencias (DS) para los cuatro modelos analizados. 


Modelo i a (240 h) a total (+) 
dp 0.637 -0.557 
Potencial 
DS 0.765 1.052 
dp 0.903 -1.05 
Exponencial 1 
DS 0.813 1.242 
dp 0.149 -0.104 
Exponencial 2 
DS 0.819 0.975 
dp 0.187 0.035 
Exponencial 3 
DS 0.758 0.912 


(*) Considera la totalidad de datos con que se cuenta; a,, coeficiente de agotamiento; í, indicador estadístico. 
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los malos resultados obtenidos tanto en las 
diferencias promedio como en la desviación 
estándar de la diferencias y en la totalidad de 
la validación estadística. 

Con el fin de profundizar aún más en el 
análisis de los modelos, se procedió a hacer una 
comparación con los promedios reales de las 
crecidas. Para ello se realizó un cociente entre 
el error estándar de estimación y los caudales 
promedio reales para las 25 crecidas estimadas, 
situación que representa la proporción del error 
frente a los caudales reales. Es decir, los valores 
cercanos a O representan la menor proporción, 
siendo estos valores los que denotan el mejor 
ajuste. Esto se presenta a continuación en el 
cuadro 4. 

Como se aprecia en el cuadro 4, fue el 
modelo exponencial 3 el que obtuvo los mejores 
resultados, es decir, el menor cociente entre el 
EEE y el promedio de los caudales observados. 

Para el caso de los modelos restantes, 
el modelo exponencial 1 y 2 no obtuvieron 
resultados satisfactorios, siendo el modelo 
de peor calidad el modelo exponencial 1, 
corroborando nuevamente, no ser un modelo 
adecuado para la cuenca en estudio. 

En relación con el lapso que dio los mejores 
resultados, en general fue el correspondiente al 
total de datos con que se cuenta. En contraparte, 
no existen evidencias notorias sobre si algún 
lapso en particular posee malos resultados 
según los estadísticos antes señalados. Sin 
embargo, al observar el test de Bland y Altman, 
es posible visualizar que el lapso con los 
mejores resultados fue el correspondiente a las 
240 horas (10 días). 


Conclusiones 


La importancia de una correcta modelación 
de este tipo de caudales, en especial en zonas 
mediterráneas, áridas y semiáridas, hace a 
este tipo de estudios de vital importancia 
en la planificación en escenarios de eventos 
extremos, como lo es la escasez hídrica, pues 
permite calcular el volumen disponible luego 
de una crecida y regular el funcionamiento de 
un embalse, entre otras acciones. 

De los resultados presentados se puede 
concluir que el modelo exponencial 3 obtuvo 
los mejores resultados para la estimación de 
caudales recesivos en la cuenca del estero 
Upeo, en Upeo, Chile. 

En la determinación del mejor lapso para el 
cálculo del parámetro a, se puede desprender 
que, en general, correspondió al total de 
datos de cada crecida, lo que era previsible, 
obteniendo sustancialmente los mejores 
resultados. Sin embargo, para el test Bland y 
Altman, el mejor lapso a utilizar resultó ser el 
de 240 horas. Dada la reconocida potencia del 
test ACBA, se concluye que el mejor periodo de 
estimación es el correspondiente a 240 horas. 

Finalmente, el modelo propuesto debiese 
responder al comportamiento de cauces y 
cuencas similares, y que en el caso de Chile 
son muchas las que comparten similares 


características geomorfológicas. 
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